Tuning effective interactions close to the critical point in colloidal suspensions 
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We report a numerical investigation of two colloids immersed in a critical solvent, with the aim of 
quantifying the effective colloid-colloid interaction potential. By turning on an attraction between 
the colloid and the solvent particles we follow the evolution from the case in which the solvent 
density close to the colloids changes from values smaller than the bulk to values larger than the 
bulk. We thus effectively implement the so-called (+,+) and (— , — ) boundary conditions defined 
in field theoretical approaches focused on the description of critical Casimir forces. We find that 
the effective potential at large distances decays exponentially, with a characteristic decay length 
compatible with the bulk critical correlation length, in full agreement with theoretical predictions. 
We also investigate the case of (+, — ) boundary condition, where the effective potential becomes 
repulsive. Our study provides a guidance for a design of the interaction potential which can be 
exploited to control the stability of colloidal systems. 



I. INTRODUCTION 

Interactions between colloidal particles in dispersions 
depend, beside the colloid coordinates, on the degrees of 
freedom of the solvent and of the co-solutes. When a 
clear separation in time and space scales between solvent 
molecules and colloidal particles exists, it is possible to 
integrate out the solvent (and co-solute) degrees of free- 
dom and derive an effective potential that describes the 
interaction between colloids [1 j. A well known case is the 
depletion interaction derived long time ago by Asakura 
and Oosawa [2] which have shed light on the role of en- 
tropic forces on the phase diagrams of colloidal suspen- 
sions. Depletion interactions arise whenever small ad- 
ditives, such as polymers or surfactants, are added to 
colloidal dispersions; when two colloids are closer than 
the typical size of the co-solute, the latter is excluded 
by a "depletion" region in between the two colloids. As 
a result, a pressure gradient originates, giving rise to a 
short-range entropy-driven attractive effective potential 
- the depletion potential. Today we know that depletion 
forces can be strong enough to induce colloidal phase sep- 
aration [3]-[5] , or to enhance the stability of the crystalline 
phase, as in the case of proteins HUE]- 

An interesting effective potential is the one arising 
when two colloids are dispersed in a critical fluid. Close 
to the critical point, thermal fluctuations of the order pa- 
rameter are correlated over length-scales that are much 
larger than the solvent size and their properties become 
independent of the microscopic behavior of the system, 
i.e. they are universal. The confinement of the order- 
parameter fluctuations in between the two colloids give 
rise to long-range effective forces, named critical Casimir 
forces [8]. 

In the last decade direct experimental evidence of such 
forces has been provided [9] . Theoretical investigation of 
critical Casimir forces, which started with the pioneering 
work of Fisher and de Gennes [10] back in 1978, has pro- 
vided accurate predictions for the radial dependence of 
the effective potential and the associated scaling proper- 
ties. Such studies, capitalizing on the universality aspect 



of the problem, have been mostly focusing on Ising model 
calculations [TTTH4] and, more recently, on a field-theory 
approach for a classical binary mixtures [T5HT9] . The 
main result of these studies is the explicit derivation of 
the critical Casimir potential for a solvent confined be- 
tween the surfaces of two large colloids: 

mz) = ^G(z/0. (1) 

where <j c is the colloid diameter and z is the distance 
between the surfaces of the two colloids. The scaling 
function Q(z/£) depends on the distance from the critical 
point (via the correlation length £), on the bulk univer- 
sality class of the solvent and on the boundary conditions 
(BC) imposed by the colloidal surface properties. It has 
been demonstrated [15] that for z » £ >> a s (the lat- 
ter being a measure of the solvent size), the behavior of 
the scaling function in Eq. ([TJ is 

e(*/0(±,±)(*/£ » i) = TA (±)±) (z/0e-<*/« (2) 

where the (=L,±) signs are related to different BCs, i.e., 
to the different absorption preferences of the confining 
surfaces with respect to the solvent: (+,+) or (— ,— ) 
corresponds to symmetric BC and (+, — ) (or equiva- 
lently (—,+)) to asymmetric BC. Experimental results 
have shown that it is possible to generate repulsive and 
attractive critical Casimir forces between a colloid and 
a substrate by modifying the BC of the substrate [9]. 
Moreover it has been shown that is possible to continu- 
ously tune the force from attractive to repulsive by pro- 
ducing a gradient in the physico-chemical properties of 
the substrate [20]. This can be exploited for inducing 
orientationl-dependent effective forces in colloids close to 
chemically patterned surfaces [2lJ [22] . Critical Casimir 
forces have also been observed when the critical behavior 
of the host medium refers not to the solvent itself but to 
the presence of small interacting co-solutes added in so- 
lution. If the inter-cosolutes interaction becomes strong 
enough, a phase separation (the analog of the gas-liquid) 
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takes place in which the co-solutes partition themselves 
in two phases of different concentration. Close to the 
corresponding critical point, the critical fluctuations in 
co-solute concentration generate critical Casimir forces. 
Such mechanism has been exploited in the experimental 
study of Buzzaccaro and coworkers [18 j, in which PMMA 
colloidal particles are dispersed in an interacting micellar 
solution. In the same work, the authors presented an in- 
teresting connection between depletion forces generated 
by the presence of the micelles far from the critical point 
and Casimir forces generated by the critical fluctuations 
close to the micellar critical point. 

Most of the theoretical and numerical investigations of 
Critical Casimir forces have been based on studies of lat- 
tice models, exploiting the universality properties of the 
phenomenon. In a recent study [23], we have reported a 
numerical evaluation of the effective interaction potential 
between two spherical hard-sphere colloidal particles im- 
mersed in a critical depletant, modeled as short-ranged 
square- well attractive potential, with the aim of exploring 
how the interaction potential changes from the depletion 
shape occurring far from the critical point, to the uni- 
versal shape induced by critical fluctuations close to the 
depletant critical point. We have confirmed the critical 
nature of the effective potential close to the critical point 
by showing that the quantity £ entering in the effective 
potentials (see Eq. [2| is consistent with the bulk corre- 
lation length of the critical co-solute extracted from the 
static structure factors. Interestingly enough, we have 
found that the strength of the effective potential between 
the two hard-sphere colloids (when the colloid-depletant 
interaction is only controlled by excluded volume, and 
hence in the (— , — ) BC class) is sufficient to drive bulk 
phase separation of the colloidal solution well before the 
critical region is approached. 

In this article we explore the changes in the effective 
potential when the colloid-cosolute interaction is continu- 
ously modified from hard-core repulsion to strong attrac- 
tion, continuously moving from the (— , — ) to the (+, +) 
BC. We also explore the interesting case in which (+, — ) 
BC are present, i.e. the case of two colloids interacting in 
different ways with the cosolutes, again interpolating be- 
tween the (— , — ) and the (+, — ) limits on changing (this 
time only for one of the two colloids) the colloid-cosolute 
interaction strength. We calculate numerically the effec- 
tive potential both at a high T, where critical phenomena 
are absent as well as close to the co-solute critical point. 



II. MODELS AND METHODS 

The total colloid-colloid interaction potential <pcc re- 
sults from the sum of the bare colloid-colloid interaction 
Vcc an d of the effective potential V e ff arising from the 
integration of the solvent and cosolute degrees of free- 
dom, 



We model Vcc as a excluded volume interaction be- 
tween two colloids of size a c 



V C c(ri 2 ) 



oo, ri2 < cr c 
ri2 > cr c 



(4) 



Inspired by the work of Buzzaccaro and coworkers [18] , 
we model the critical medium as a fluid of interacting co- 
solutes dispersed in an implicit solvent. The co-solutes, 
of size <j s interact via a pairwise square-well potential 
(SW) 



V S s (Uj) 



OO, Ti 

o n 



< Tij < (1 + S)(T S 

> (1 + S)a s 



(5) 



^cc{r 12 ) = V C c(r 12 ) + V ef f(r 12 ) 



(3) 



where e s controls the strength of the interaction and S the 
relative (respect to a s ) width of the well. a s and e s are 
chosen as unit of length and energy. The temperature T 
is measured in units of e s . The co-solute is characterized 
by a gas-liquid critical point located at (T c = 0.478, <p c = 
0.25) [24 , where (j) c = (7r/6)p c <7g is the critical packing 
fraction of the co-solute and p c is its number density. For 
this model, it has been shown [23] that the correlation 
length and the susceptibility extracted from the static 
structure factors close to the critical point diverge with 
a power-law with the respective Ising critical exponents. 
The size ratio q = cr s /a c between the co-solute and the 
colloids is fixed at q = 0.1. 

To evaluate the effective potential, we perform Monte 
Carlo simulations of two colloids in a fluid of co-solute 
particles, in the canonical ensemble at fixed T and p s in 
a parallelepiped shape. The two colloids are constrained 
to move only along the x— axis, sampling only a limited 
range of distances. Several overlapping relative distances 
windows are simulated, evaluating for each window P(r), 
the probability of observing the two colloids at relative 
distance r. Splicing together the P(r) evaluated in differ- 
ent windows provides an effective (and parallel) way for 
evaluating the entire P(r). The logarithm of P(r) is by 
definition the effective potential (apart from an overall 
constant which is fixed imposing V e ff(oo) = 0). To min- 
imize finite-size effects at the temperatures investigated, 
the dimensions of the box (L x ;L y ;L z ) are chosen in such 
a way that the surface-to-surface colloidal distance eval- 
uated via the boundary conditions is more than twice the 
distance over which the effective potential goes to zero. 
Moreover, along all directions, the solvent density profile 
reaches a constant value on approaching the box bound- 
aries. Close to the critical point and along the critical iso- 
chore, the size of the box is [L x = 52; L y = 26; L z = 26], 
requiring 16000 co-solute particles. The bulk density is 
estimated a posteriori by calculating the local density far 
from the two colloids. 

We also evaluate the co-solute density profile for differ- 
ent BC by selecting a volume centered along the x— axis 
of transversal section equal to a 2 s and we average the local 
density with a mesh of the order of 0.1<r s . 
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In this paper, when discussing the critical behavior, 
we will show results for effective potentials and density 
profiles evaluated at the critical packing fraction and at 
the reduced temperature T/T c = 1.0251, corresponding 
to a critical correlation length £ = 2.ba s [23] . 



III. RESULTS 

A. (--) BC 

In our previous work [23] we have discussed the evo- 
lution of the effective potential when the interaction be- 
tween the colloids and the solvent Vcs are treated as 
hard-spheres. In such condition, close to the colloid sur- 
face, the density of the solution is smaller than the aver- 
age, effectively generating a (— , — ) BC. We have shown 
that close to the critical point, V e ff is long-ranged, sig- 
naling the onset of criticality. Its radial dependence is 
well described by the exponential decay of Eq. [T] with 
the same characteristic length £ of the bulk critical fluc- 
tuations. 

To get more insight into the mechanism which drives 
attraction between colloids close to the critical point in 
the (— , — ) BC case, we show in Fig.[l]the co-solute den- 
sity profile along the horizontal x-axis for two different 
relative colloid-colloid distances. Notice that the two col- 
loids are always located symmetrically with respect to the 
origin. Thus the density profile is symmetric with respect 
to x = 0. For this reason in Fig. [T] we show only the pos- 
itive x region. In between the two colloids, the density 
is significantly lower than the bulk density, a typical de- 
pletion effect. Outside, the density relaxes toward the 
bulk value with an exponential decay, again controlled 
by the bulk critical correlation length. It is interesting 
to discuss the physical origin of the net attractive force 
between the two colloids. Since the co-solute colloid in- 
teraction is modeled via an hard-sphere potential, it is 
possible to prove that the effective force results from the 
mismatch in the contact density along the two sides of 
each colloid. Fig. [I] shows that indeed, at contact, the 
density outside is slightly larger than inside. The figure 
also shows that the mismatch decreases on increasing the 
relative distance between the two colloids in parallel with 
the decrease of the effective force. 



B. From the (-, -) to the (+, +) BC 

To drive the transition from (— , — ) toward (+, +) BC 
we tune the strength of a short-range attraction in the 
colloids-cosolute interaction Vcs- The attraction gives 
rise to an enhanced accumulation close to the surface 
of the colloids, at first compensating and then inverting 
the depletion characteristic effect of the excluded volume 
interaction [25] [26] . 

We model the attraction via a pairwise SW potential 
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FIG. 1. Normalized co-solute density profiles with re- 
spect to the bulk density po along the horizontal x-axis at 
T/T c = 1.0251 and at the critical density C = 0.2501. Co- 
solute particles interact with a short-range SW attraction 
while colloid-cosolute interaction is modeled with a HS repul- 
sion. The two profiles are evaluated for two different colloid- 
colloid distances, namely for r = <j c + 3cr s and r = a c + 4cr s , 
with cr c = 10cr s . The box is centered at x = and the two col- 
loids are centered at ±r/2 along the x-axis. Being the profiles 
symmetric with respect to x = only positive x are shown 
here. Note that each colloid excludes the solvent in a region 
equal to <j c + cr s . The dashed lines highlight the mismatch in 
the density at contact between the inner and outer side of the 
colloid. The inset shows that the outer part of the density 
profile is well fitted by an exponential p = po + Aexp(— r/£), 
where £ has been fixed to £ = 2.5cr s , i.e. to the bulk critical 
correlation length value for T/T c — 1.0251 at <fi c . In this case 
the fit parameter A < 0. 



{oo, r < a cs 
-e, a cs <r < a cs + Aa s (6) 
r > a cs + Aa s 

where <j cs = (<j s +<j c )/2 and A = 0.35. The width A<j s 
has been chosen as a compromise between limiting the 
colloid-cosolute interaction to the nearest-neighbor shell 
and maximizing the volume over which cosolute bind to 
the colloid. The parameter e is used as control parameter 
to drive the cross-over from the hard-sphere like behavior 
(s — » 0) to the wetting case s/kT > 1). 

We start by discussing the behavior of the effective 
potential at high T, where critical phenomena are not 
present. The evolution of the effective potential upon 
changing e is shown in Fig. [2] When e = 0, the well 
known depletion interaction potential is observed. On in- 
creasing £, the depletion attractive interaction is progres- 
sively weakened, and the potential at contact becomes re- 
pulsive. For even large e values, the colloid becomes sur- 
rounded by a persistent layer of co-solutes which extent 
the effective radius of the colloid, making it impossible 
to attract a neighboring colloid for distances closer than 
<j c + a s [27]. Under such strong coupling conditions, the 
effective potential acquires an oscillatory character, with 
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minima originating from the preferential distances allow- 
ing for an integer number of co-solute layers between the 
colloids. The configuration associated to the first mini- 
mum, called bridging, is the most energetically favorable 
since one co-solute particle is bonded with both colloids. 
The other minima are related to particular configura- 
tions in which bridging is obtained by particle chains. A 
sketch of such situation is shown in Fig. [2] The inter- 
mediate maxima occur when the co-solute particles (or 
particle chains) are not bonded to both colloids. For ex- 
ample the first maximum corresponds to a situation, also 
illustrated in Fig. [2j in which a single co-solute particle 
cannot be bonded to both colloids, since r > a s . The evo- 
lution of the effective potential reported in Fig. [2] clearly 
show how the minimum at contact progressively turns 
into a maximum, and simultaneously a new minimum 
develops at the bridging distance. Intermediate values of 
s thus provide a viable mechanism for contrasting deple- 
tion interaction and favoring colloidal stability. The case 
e = 0.5 is emblematic, since the minimum at contact has 
essentially disappeared while the new minimum at dis- 
tance 1.1<t c has not yet developed, such that the effective 
potential is never (in absolute value) significantly larger 
than ksT. For very large e values, V e ff becomes again 
sufficiently intense to drive a colloid phase separation. 
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FIG. 2. Effective potential for two large colloids inter- 
acting through a hard-core repulsion at T/T c — 2.09 and 
4> — 0c = 0.25. Colloids are in solution with co-solute par- 
ticles interacting (among themselves) via a SW attraction 
of width 5a s = 0.1a s and depth e s . The two colloids in- 
teract with co-solute particles via a SW attraction of width 
Acr s = 0.35cr s and depth e. Sketches in the figures represent 
particular configurations adopted by the co-solute for differ- 
ent colloid-colloid distances. Such configurations give rise to 
maxima and minima in the potential. 

Close to the critical point, critical Casimir forces add 
to the standard effects described before. Differently from 
the pure (— , — ) case, the presence of an attractive inter- 
action between the colloid and the solvent brings in two 
new features that are distinguishable for medium and 
large values of e: a repulsive contribution to the effective 
potential at short distances and the presence of oscilla- 
tions induced by the granularity of the cosolutes. The 
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FIG. 3. Effective potential of two large colloids close to the 
critical point of the co-solute (T/T c = 1.0251 and cj) = <frc = 
0.25). The two colloids interact with co-solute particles via 
a SW attraction of width Aa s = 0.35cr s and depth e. The 
two colloids have therefore (+,+) BC. The increase of the 
attraction strength give rise to large oscillations associated to 
bridging effects as discussed in Fig. [2] 

theory on critical Casimir forces predicts that at dis- 
tances larger than £ the effective potential develops an 
attractive exponential tail, so that the effective potential 
behaves the same for the (— , — ) BC and the (+, +) BC 
cases. Fig. [3] shows how the effective potential evolves on 
increasing the value of e. As expected the contact value 
becomes repulsive but the long tail behavior remains al- 
ways attractive and has the same exponential character 
of the e = case [23 . Even for large values of e, where 
the same oscillations characterizing the high T effective 
potentials modulate the shape, the overall behavior can 
be represented by an exponential function. The density 
profile for the (+, +) case is shown in Fig. [4] Close to the 
colloidal surfaces a significant layering of the cosolutes is 
observed, consistent with the presence of oscillations in 
V e ff(r). On the external sides, beyond the oscillations, 
the density decays again with an exponential shape con- 
trolled once more by £. In comparison with the (— , — ) 
case, the resulting interaction between the two particles 
is more difficult to visualize, since it arises from the com- 
petition between the density at contact (in which the sol- 
vent pushes the colloid) and the density the well bound- 
ary (where the solvent attracts the colloid). Indeed, in 
the case of a colloid of diameter sigmac interacting via 
square- well attraction with solvent particles, the pressure 
originates from the two points of discontinuity of the po- 
tential, the colloid-cosolute hard-wall distance (HWD) 
and the square-well distance (SWD) a c s + Aa s accord- 
ing to the expression [28] : 

P pCX 2 r 3 , +x 
—— = ~Kp[a cs g((7^ s ) 
P 3 

- (a cs + Aa s ) 3 g(cr cs + Aaj)(l - e^)}. (7) 
where g(cr cs ) is the pair correlation function evaluated at 
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the HWD and g(a cs + Aa s ) at the SWD. From Eq.[F]one 
can notice that the density at the hard-wall generates a 
positive contribution to the pressure, while the density 
at the well boundary provides a negative contribution. 
Figure [4] shows that indeed the density at the SWD in- 
ner side (left side of Fig. [4| is significantly higher than 
the density at the SWD outer side (right side of Fig. [4| , 
and it is responsible for the resulting net attractive force 
between the two colloids. Indeed, the excluded volume 
contribution in this case would tend to separate the two 
colloids, being the contact density outside slightly larger 
than the contact density inside. 
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FIG. 4. Normalized co-solute density profile with re- 
spect to the bulk density po along the horizontal x— axis at 
T/T c = 1.0251 and at the critical density C = 0.2501. Co- 
solute particles interact with a short-range SW attraction of 
width 5a s = 0.1a s and depth e s while colloid-cosolute interac- 
tion is modeled with a SW interaction of width Acr s = 0.35a s 
and depth e = le s . The box is centered at x = and the 
two colloids are centered at dbr/2 along the x-axis. Being the 
profile symmetric with respect to x = only positive x are 
shown here. Note that each colloid excludes the solvent in a 
region equal to <j c + cr s . The dashed segments highlight the 
mismatch in the density at contact (HWD) and at the well 
boundary (SWD) between the inner and outer side of the col- 
loid. The profile is associated to the effective potential in Fig. 
3 (squares). As for Fig. 1, the criticality of the co-solute can 
be caught from the exponential decay of the density outside 
the colloids towards p s (x)/po = 1. The exponential fit shown 
in the figure (solid line), has been performed by fixing the 
exponential decay with the value of the bulk critical correla- 
tion length for the temperature and density investigated, i.e. 
£ = 2.5a s . 



C. From the (- 



to the (+,-) BC 



We now discuss the situation in which one of the 
two colloids (Ci) interact with the co-solute through a 
square- well attraction (Vc 1 s(^u)) while the other col- 
loid (C 2 ) experience only excluded volume interactions 
(Vc 2 s(r2i)). More precisely 



V Cl s(ru) 




(8) 



and 



t 7- / \ J 00, r 2j < (J cs 
10 r 2j > cr cs . 



(9) 



Also in this case we start investigating a T signifi- 
cantly larger than the critical one for different interaction 
strengths e. The results are shown in Fig. [5] 
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FIG. 5. Effective potentials for two large colloids interacting 
through a hard-core repulsion at T/T c = 2.09 and cj) — <j) c — 
0.25. Co-solute particles, in solution with the two colloids, 
interact through a SW attraction of width 5 = 0.1a s and 
depth e 8 . One of the two colloids (Ci) interact with co-solute 
particles via a SW attraction of width A = 0.35cr s and depth 
e. The other colloid (C2) has instead a HS repulsion with the 
co-solute. 

The trend is similar to the one observed going from 
( — , — ) to the (+, +) case; (i) the contact value of the 
effective potential grows continuously on increasing e, 
progressively offsetting the original e = depletion in- 
teraction. Still, the strength of the repulsion is signifi- 
cantly smaller than the one observed in the (+, +) case, 
(ii) only weak oscillations (with amplitude smaller than 
fcfiT) characterize the radial dependence of the effective 
potential, signaling the absence of a strong layering of the 
cosolute between the colloids. We notice that similar fea- 
tures have also been observed in the case of non-additive 
HS mixtures 29) [30]. 

Close to the critical point, the theory predicts that, on 
increasing e, V e ff changes from an attractive to a repul- 
sive exponential decay, on changing the boundary condi- 
tions from (— ,— ) to (—,+). The effective potentials for 
different e close to the critical point are plotted in Fig. [6] 
While for r < 1.3<j c the echo of the layering effects is 
still visible, for larger r values the decay of all curves 
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FIG. 6. Evolution of the effective potential of two colloids 
close to the critical point of the co-solute (T/T c = 1.0251 and 
cj) — cj) c — 0.25) on increasing e. One of the two colloids (C\) 
interact with co-solute particles via a second SW attraction 
of width A = 0.35cr s and depth s. The other colloid (C2) 
interacts via a HS repulsion with the co-solute. This setting 
provides a realization of (+, — ) BC. 



is compatible with the same exponential function decay, 
again supporting the identification of the interaction po- 
tential in this spatial region as arising from the universal 
behavior imposed by the critical fluctuations. As in the 
case of (— , — ) BC, for (+, — ) BC the density profile pro- 
vides useful informations on the mechanism that gives 
rise to the repulsion. Fig. [7] shows the co-solute particles 
density profile along the x— axis. The different colloid- 
cosolute interaction result in a different density profile 
around the two colloids. Around C2 (— BC), the solvent 
density is lower than the average, while the opposite be- 
havior is observed for C\ (+ BC). The repulsive force 
on C2 originates from the slight mismatch of the contact 
density (larger inside than outside), while the repulsive 
force on C\ originates from larger density at the SWD on 
the outside, compared to the inside. Contrary to what 
observed in the (— ,— ) to (+, +) BC case, here tuning 
the attraction strength allows us to modify the sign of 
the effective force [9]. In fact we observe that for small 
e the effect of the colloid-cosolute attraction is not suf- 
ficiently strong to substantially change the shape of the 
effective potential, which remains completely attractive. 
For higher values of e the potential turns into a com- 
pletely repulsive one and no sign of oscillations driven by 
co-solute structures is visible. It is interesting to note 
that under these conditions it is in principle possible to 
tune finely e in order to obtain a flat co-solute density 
profile (apart from the layering at contact). This corre- 
sponds to impose Dirichlet BC [12] and from a conceptual 
point of view to identify a sort of 6 condition, in analogy 
with polymer solutions pQ , where the effective interaction 
potential is close to zero. 




■ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

-20 -10 10 20 

x 



FIG. 7. Normalized co- solute density profile with re- 
spect to the bulk density po along the horizontal x— axis at 
T/T c = 1.0251 and at the critical density <j) c — 0.25. Co- 
solute particles interact with a short-range SW attraction of 
width 5 = 0.1a s and depth e a . One of the two colloids (Ci) in- 
teracts with co-solute particles via a SW attraction of width 
A — 0.35cr s and depth e — ls s . The other colloid(C2) in- 
teracts via a HS repulsion with the co-solute. The profile 
is associated to the effective potential in Fig. [6] (squares). 
The different colloid-cosolute interaction results in a different 
density profile around the two colloids. Around C2 (colloid 
on the right), the solvent density is lower than the average, 
while the opposite behavior is observed on C\ (colloid on the 
left). The repulsive force Ci colloid originates from the slight 
mismatch of the contact density (larger inside than outside), 
while the repulsive force on C2 originates from larger density 
at the square- well location on the outside, compared to the 
inside. The exponential fits (solid lines) have been done by 
setting the exponential decay to the value of the bulk critical 
correlation length £ = 2.5cr s . 



D. From the (+, +) to the (+, -) BC 

According to theoretical predictions [15] the effective 
potential is expected to change in the transition from 
(+, +) to (+, — ) exactly as in the case from ( — , — ) to 
(+,—). These predictions refer to the scaling region of 
the potential. In the present numerical study, we also 
access the short distances, where the effective potential 
probes the non-universal aspects of the solvent-colloid in- 
teraction. To highlight the difference at short distances 
between the two cases, Fig. [8] shows the evolution of the 
effective potential when the attraction between only one 
of the two colloids and the co-solute is progressively re- 
duced. Differently from the results of Fig. [6] in this case 
only the long distance behavior changes sign, while the 
short distance part of the potential remains always repul- 
sive. Finally in Fig. [9] we plot cases differing in their BCs, 
but all at the same temperature and co-solute critical 
packing fraction. In all cases, the long-distance behav- 
ior of V e ff can be described, as expected theoretically, by 
an exponential decay, with a correlation length consistent 
with the bulk critical correlation length £. 
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FIG. 8. Evolution of the effective potential for two colloids 
close to the critical point of the co-solute (T/T c = 1.0251 and 
cj) — cj) c — 0.25). One of the two colloids (Ci) interacts with 
co-solute particles via a SW attraction of width A = 0.35cr s 
and depth s\ — 1.0£ s . The other colloid (C2) interacts with 
the co-solute through a SW attraction as well with the same 
width A but depth £2. Starting from £2 = £1, i.e. (+, +) BC, 
the interaction strength £2 is decreased down to the limit case 
of £2 = corresponding to (+, — ) BC. Contrary to the case 
illustrated in Fig. [6] here only the long distance tail can be 
tuned from attractive to repulsive, while the short range part 
of the potential remains repulsive. 




FIG. 9. Effective potentials for different BC, evaluated at the 
co-solute critical packing fraction and at T/T c — 1.0251. In 
the critical region, independently on the BC, the long range 
part of the effective potential decays exponentially. According 
to eq. |l]) it is possible to extract the characteristic length that 
controls the decay of the exponential. To extract this value 
the three curves have been interpolated in the same range with 
eq. |l]) starting from r = 1.3cr c . For the three cases we find 
£(-,-) = 2.24cr s , £(+,+) = 2.49cr s , and £(+,-) = 2.20cr s which 
is consistent the bulk correlation length £ = 2.5cr s of the co- 
solute for this T/T c 23 . The non universal amplitudes found 
from the exponential fits are = —0.63 = —0.45 

0.48. 



IV. CONCLUSIONS 

In this article we have discussed how the effective po- 
tential, resulting from confining critical co-solute parti- 
cles between the surfaces of two large colloids, depends on 
the colloid-cosolute interaction. In a previous study [23] 
we have already investigated the case of (— , — ) BC, i.e. 
when only colloid-cosolute hard-core repulsion is present. 
In such case we have shown that the resulting potential 
for T/T c — >• 1 at the co-solute critical packing fraction is 
monotonic, attractive and long range. The co-solute den- 
sity profile along the x— axis of the simulation box shows 
that for (— , — ) BC the density close to colloids is lower 
than the bulk density and it relaxes exponentially to the 
bulk value far from the two colloids. Upon switching on 
the colloid-cosolute attraction e , depletion effects are 
progressively weakened and the co-solute density close to 
the colloid varies from values smaller than the bulk to 
values larger than the bulk. Correspondingly, the con- 
tact value of the potential varies from negative to posi- 
tive while the long distance part of V e ff (r), dominated 
by critical fluctuations, remains always attractive as the- 
oretically predicted in the (+, +) BC case. The increase 
in s has instead a profound effect on the non-universal 
short-distance part of the potential, which progressively 
develops large oscillations. Such oscillations are related 
to specific geometries at characteristic lengths associated 
to integer number of co-solutes between the two colloids. 
At these specific distances, the local energy is minimized. 
For intermediate values of e the colloid-cosolute attrac- 



tion compensates the depletion interaction, effectively re- 
ducing the short-distances interaction potential. We have 
also addressed the case of two colloids interacting differ- 
ently with the co-solutes, a realization of the so-called 
(+, — ) BC case. In this case, the local density close to the 
two colloids is respectively lower and higher than the bulk 
value and the critical Casimir forces are expected to be 
repulsive. We have shown that indeed, when the asym- 
metry in the interaction potential is sufficiently intense to 
drive different BC, the potential becomes repulsive at all 
length scales. We have shown that, independently on the 
BCs, the critical long range part of the effective potential 
is always described by an exponential whose decay is con- 
trolled by the critical correlation length £ corresponding 
to the thermal correlation length of the bulk co-solute 
close to the critical point, in full agreement with theoret- 
ical predictions [15]. We have investigated the behavior 
of the density profile in all examined BC cases. Building 
on the fact that all interactions are modeled as excluded 
volume or as square- well attraction, the net pressure over 
the colloids can be estimated simply by the density of co- 
solute at contact and at the square-well distance (in the 
case of attraction between the colloid and the cosolute). 
The density profiles confirm that in the (— , — ) case, the 
contact density inside is smaller than outside, determin- 
ing the net attraction. In the (+, +) case a different 
mechanism for attraction is observed: the leading con- 
tribution arises from the mismatch of the density at the 
well distance, larger inside than outside, determining a 
net attraction. Finally, in the (+, — ) case, repulsion is 
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driven by two different mechanisms for the two colloids. 
The hard-sphere colloid is pushed out by the larger con- 
tact density inside. The attractive colloid is dragged out 
by a larger density at the square well distance. 

The possibility of varying the potential from repulsive 
to attractive and to finely control its shape and intensity 
by tuning the BC, provides a guidance for controlling 
equilibrium properties of colloids dispersed in precriti- 
cal suspensions. This is important for future applica- 
tions, for instance in the case in which the confining sur- 
faces are chemically patterned colloidal particles, such as 
patchy [31] or Janus particles [21] [32] . In this case, the 
introduction of a geometrical constraint in the colloid- 



cosolute interaction would give rise to a torque that can 
be used to control the orientation of colloids, as shown 
for chemically patterned substrates [22], and that could 
give rise to new unexplored phases. More in general, our 
results provides useful informations for designing inter- 
action potentials which can be exploited to control the 
stability of colloidal systems. 
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